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Abstract 

The coupled cluster or exp S form of the eigenvalue problem for lattice Hamiltonian 
QCD (without quarks) is investigated. 

A new construction prescription is given for the calculation of the relevant coupled 
cluster matrix elements with respect to an orthogonal and independent loop space ba- 
sis. The method avoids the explicit introduction of gauge group coupling coefficients by 
mapping the eigenvalue problem onto a suitable set of character functions, which allows 
a simplified procedure. 

Using appropriate group theoretical methods, we show that it is possible to set up the 
eigenvalue problem for eigenstates having arbitrary lattice momentum and lattice angular 
momentum. 



1 Introduction and Overview 



The investigation of the eigenvalue problem for the lattice QCD Hamiltonian is considered to 
be an alternative to standard Lagrangian lattice Monte Carlo QCD, possibly giving new insight 
into the structure of such non-abelian gauge theories. 

For pure SU(3) Yang-Mills theory |L| (without fermions) in particular, many attempts have 
been made to attack the corresponding Kogut-Susskind Hamiltonian problem: for instance, 
there exist the strong coupling expansion [fj], the exp(-tH) method^] or variational techniques 0. 
Up to now, none of these approaches could obtain results for excited states (e.g. glueball masses) 
comparable in control and accuracy to those within the Euclidean Monte Carlo method (there 
has been, however, some progress for ground states using the Greens function Monte Carlo 
method 0). 

This also holds for the coupled cluster (exp S) method which attracted special attention 
in recent years |], |7], ^, [|. (Some encouraging results within this framework were obtained 
recently [K|.) Here the basic idea is to incorporate manifestly the correct volume dependencies 
of observables by writing the ground state in the form ipo = e s and putting ip = Fipo for excited 
states. The " Schroedinger" equation for the functions S and F can be formulated rigorously || 
and it is tempting to define approximations by a suitable truncation of a loop space expansion 
of these quantities |], 0, |10| . 

It is the purpose of this paper to further elucidate the structure of this coupled cluster 
method with the hope that the resulting insights may lead to improved calculations of the 
QCD spectrum. 

We will concentrate our considerations on the treatment of the Kogut-Susskind Hamiltonian 
as the lattice regularization of an SU(n) Yang Mills theory. A discussion of the full QCD and 
its treatment within a quenched approximation is possible, but this will be deferred to a future 
publication. 

We now give an outline of our paper which summarizes at the same time our methods and 
our results. 

Our basic tools will be group theoretical methods which will be introduced in section 2. 
The group of the link variables, the local lattice gauge group and the lattice Euclidean group 
will play a role. 

As discussed in section 3, projection operators on representations of the lattice Euclidean 
group with given lattice momentum and lattice angular momentum allow one to introduce the 
notion of an intrinsic wavefunction related to the ground state function S and to the "excitation 
operator" F. 

This structure has been used in Refs. [||, |7|, |l(J for the trivial representation; here we provide 
a systematic framework for general representations of the lattice Euclidean group. 

The solution of the eigenvalue problem for the Kogut-Susskind Hamiltonian is then reduced 
to the determination of the intrinsic eigenfunctions. 

For this purpose, a basis of suitable wavefunctions is needed which may be used for an 
expansion and which allows a computation of the relevant coupled cluster matrix elements. 
Within the Kogut-Susskind theory these have to be functions of the link variables which are 
invariant under the action of the local lattice gauge group. 

The problem of setting up such a basis in an effective way is addressed in section. 4. There 
exist two strategies for the construction of such basis systems: 

1) Choose first a basis for the functions of the individual link variables given by the standard 
D-functions. General polynomials of these functions with different link variables, combined with 
suitable SU(n) coupling coefficients, form then the desired basis for the intrinsic hadron (or 
vacuum) wave functions. We call this set of functions the D-loop basis. 



Details of this construction have been worked out in Ref.jll]. An application is the "exact 
linked cluster expansion" discussed in Ref.||12|. 

This method is limited by the necessity to handle an increasing number of SU (n) coupling 
coefficients. 

A clear merit of the procedure is that it provides an independent, orthogonal and (in the 
limit of increasing polynomial degree) complete basis of physical states. 

2) An alternative system of physical states is provided by the set of character functions cor- 
responding to an expansion in terms of suitable group characters. This approach was used in 
the recent calculations within the coupled cluster method|| [7|. The obvious advantage here is 
that each term is manifestly locally gauge invariant, and no coupling with SU(n) Clebsch Gor- 
don coefficients is needed. The problem, however, is that the emerging system of wavefunctions 
is in general non-orthogonal and overcomplete. 

In Refs.[|G], 0, H the disease of having linear dependencies was cured with the help of a 
special form of the Cayley-Hamilton relation for SU(n) matrices. This method, however, does 
not appear to be very systematic, and only calculations with wavefunctions generated from up 
to fourth order plaquette polynomials have been possible up to now. 

In section 5. we will introduce a new procedure for working with the orthogonal and 
independent D-loop basis which combines the above two alternatives by constructing a suitable 
mapping of the character functions on the D-loop basis avoiding, however, the explicit handling 
of SU(n) coupling coefficients. In this framework, the Cayley-Hamilton relationship in its 
general form is mainly used for systematically computing certain norm relations. 

Our procedure relies essentially on the following observations: 

1) The characteristic coupled cluster matrix elements emerge as a byproduct when the (non- 
orthogonal and overcomplete) character functions are set up systematically by an iteration 
procedure. 

2) The D-loop functions can - up to a normalization factor - be uniquely characterized by 
the eigenvalue pattern of a certain set of commuting Casimir operators. 

3) The matrix elements of these Casimir operators are computable within the character 
functions by the same methods which were used to set up these functions. 

Diagonalizing all necessary Casimir operators in the space of character functions yields then 
the mapping on the D-loop states. 

This solves the problem of linear dependencies among these functions by using the eigen- 
value patterns of the Casimir operators and by computing the relative norms of the dependent 
eigenstates with the Cayley Hamilton relation. 

The final calculations are in this way reduced to a calculation of the Kogut-Susskind eigen- 
value problem within the D-loop basis expansion. 

We hope that this will allow future numerical Hamiltonian lattice QCD calculations which 
may go to higher order than the previous attempts]^, 0, 

Also, our procedure yields a natural truncation prescription for the corresponding coupled 
cluster equations because the D-loop basis is orthogonal and unique. 

Some details of a computational strategy are described in section 6. 



2 Group Theoretical Structures 

We shall first give the definitions and notations for the SU(n) lattice Yang-Mills theory, espe- 
cially its group theoretical content. 

The general framework was given by Kogut and Susskind[|l]|. Accordingly, one has to define 



a Hilbert space Ti. given by the set of "top" wave functions depending on N link variables 



n = Mu lt ..u N )} (i) 

where the quantities U\ (/ = 1,..,N) are elements of the gauge group SU(n) and N is the 
number of oriented links in a D-dimensional lattice (D is the number of space dimensions). 

As in thermodynamics we shall work with a finite volume, i. e. with a finite lattice, for 
definiteness. However, our computational framework allows one to take an infinite volume limit 
(N — > oo) at any later stage. 

The scalar product is given by an N-fold Haar measure integral. 

The group theoretical nature of the link variables Ui gives as a natural orthogonal and 
complete basis of Ti all N-fold products of SU(n) D-functions, e.g. for SU(2) we have the 
functions 

The group of (time independent) local lattice gauge transformations is abstractly given by 

G loc = [SU(n)} M (3) 

where M is the number of sites of the lattice. 

Elements of Gi oc are written as g = g{x) where x denotes any lattice site. A unitary 
representation of Gi oc on Ti is then given by 

(p(g)^)(U 1 ,..,U N ) = y(U?,..,U 9 N ) (4) 

where the link variables are transformed like parallel transporters: 

=g(x)U l9 -\x + ee 3 ) (5) 

if the link I = (x, ej) connects the sites x and x + eej (e is the lattice spacing, ej is a positive 
unit vector in j-direction). 

The physical Hilbert space is defined by the subspace of Ti corresponding to the trivial part 
of the decomposition of the representation p, i.e. by the gauge invariant states 

n phys = {^eH\ p(g)V = * for all g G G loc } (6) 

A systematic construction of a basis of H p h ys generalizing Refs. |Tl|, 0] will be the main topic of 
this paper and is described in section 4. 

We want to impose on this basis the classification of being characterized by the irreducible 
representations of the lattice Euclidean group, which is a strict symmetry group of the lattice 
Kogut-Susskind Hamiltonian. The lattice Euclidean group is a discrete remnant of the standard 
continuum Euclidean group and is defined as follows: 

Let 

D 

R?att = {x = e n J e j I n 3 = mteger} (7) 

3=1 

be the set of lattice sites of an infinite lattice. The lattice translation group Gu is then isomor- 
phic to Rf att and given by the mapping of Rf att 

x — > x + a (8) 

for any a G R® tt - The lattice rotation group G\ r is the restriction of the group 0(D) leaving 
Rf att invariant. We call G[ r the cubic group [0, it is discrete and has 8 elements for D = 2 



and 48 elements for D = 3. The structure of this cubic group and its representations are well 
known |T3"lj. The Euclidean group Ge is then the semidirect product Ge = Gi r ® s Gu defined 



?D 



for u = (R, a) £ G E by the mapping of R\ att 

x — > ux = Rx + a (9) 

Since the mappings u may change the orientation (i.e. it may be that det(R) = — 1), the group 
Ge acts on the set of links with both orientations. We use the notation A = (I, a), a — ±1 for 
these generalized links: 

(/, 1) stands for the links with the originally chosen orientation, i.e. they have the structure 
(I, 1) = (x, ej), connecting x to x + eej (x £ Ri* att ) where ej is positive. 

(I, —1) = (x + eej, —ej) connects then x + eej to x. 

Writing A = (x, Cj) for a general link, Cj being a positive or negative lattice unit vector, the 
action of u = (R, a) £ Ge is simply given by 

A — » uX = (ux, Rcj) (10) 

This allows us to define a unitary representation T of the lattice Euclidean group Ge on 
the Kogut-Susskind wave functions as a combination of the corresponding permutation of the 
link variables and the mapping U — > U~ x if the link in question is reoriented: If \1/ depends on 
the variables U^,..,!/^ and if we put u(l a , 1) = (n a ,cr a ) (a = l,..,r, a a = ±1), then T(m)\P 
depends on the variables U ni , .., U nr and we have 

(T(«)¥)(17 B1I .., U nr ) = *(C£\ .., U£) . (11) 

As in the formal continuum limit, the operators T(u) commute with the Kogut-Susskind Hamil- 
tonian for all u £ Ge- 

We now construct projection operators on subspaces of H or Hphys corresponding to specific 
irreducible representations of Ge- 

For translations we have a "lattice momentum projection" 

n»(p)= E e- ia ^T(R=l,a) (12) 

where p £ R D is restricted to the first Brillouin zone (— n < epj < ir). 

If o% v , denote the D-functions for the irreducible representations (including inversions) of 
the cubic group G; r [|13[], a projection on "lattice angular momentum" T is given by 

n, P (r; i/,i/)= E dl u ,(R)T(R,a = 0) (13) 

By construction these definitions guarantee for any ^ £ TC the characteristic relations 

T(l,a)II,t(p)* = e ip i a m lt (p)V (14) 

A combination of both projections yields states with "good" momentum and angular mo- 
mentum in the sense that we have for 



*^ = n, t (p)n Ir (r, 



(15) 



the relations 



T(l,a)*£ = e^«^ (16) 



The basic problem of a "lattice Yang Mills theory" is then to find in H. p h ys (approximate) 
eigenfunctions of the type \E^/ of the Kogut-Susskind Hamiltonianfl| H KS = f-i? with 

H = E la E la -xV (17) 

v = Exu 



x 



9 



where g is the coupling constant and a is a colour index (a = 1, ..,n 2 — 1). Summation over 
repeated indices is always assumed; □ labels the plaquettes, and \n is given by 

Xn:=trU£U£U£U£ (18) 

when □ = (Jij <7i, . . . , Z 4 , <t 4 ). The "colour-electric field operators" E ia generate - in analogy to 
the standard momentum operator - a left multiplication of group elements in the arguments of 
the wave functions. They are quantum operators conjugate to the link operators Ui obeying 
the commutation relations 

[Eia,U v ] =^A a C/ ; , (19) 
where the SU(n) generators A a are normalized according to tr X a X b = 5 ab /2. 



3 The exp S method and intrinsic wave functions 



A motivation for the introduction of the coupled cluster or exp S method is given by the 
following considerations: 

Given a Hamiltonian H, a standard and often successful method to get the approxi- 
mate spectrum of the low lying energy states is provided by the Lanczos approach: Choose 
some trial state <fi and diagonalize H restricted to the finite-dimensional space spanned by 
(</>, H(p, H 2 <j), .., H n (f)). There are many cases where this gives reliable results if n is large enough. 



However, for our lattice Yang Mills case, this procedure is doomed to fail[|14[ because we 
have here a situation analogous to nuclear matter, for instance. In the infinite volume limit 
(N — > oo) - where we want to formulate our approach - the groundstate energy E and excitation 
energies E — E Q of H^s have the behaviour 

E oc N (20) 
E-E oc 1 

Also the groundstate wavefunction displays a characteristic "pathology" in sense that its norm 
(defined by the N-fold Haar measure integral) has an essential singularity for N — > oo. Its 
precise structure will be given below, within perturbation theory it is related to the appearance 
of disconnected diagrams. 

It has been known for a long time that this difficulty is cured by rewriting the eigenvalue 
problem within the exp S framework (see Ref.|15j for the standard many-body theory and 
Ref.[]T6] for the Kogut-Susskind theory). 



For our case, the method consists of introducing the ansatze 

* {U)=e*pS{U) (21) 
(U = (U\, .., Un)) for the ground state and 

V(U) = F(U)expS(U) (22) 

for excited states. 

The mentioned "pathology" of the ground state consists then in the fact that we have the 
norm relation \S\ 2 oc iV for the function S(U) appearing in the exponent with respect to \l/o! 

The validity of these volume dependencies is related to a characteristic linked cluster struc- 
ture of S(U) and F(U), which follows from rewriting the Schrodinger equation in terms of these 
functions, resulting in the non-linear equation 

+ SpSp -xV = E (23) 

for S and in the linear equation 

F^ + 2S^ = (E-E )F . (24) 

for the excitation operator F. 
Here, we use the abbreviation 

M=M (25) 

and the notation 

U=[E la J] , fw=[E la ,[Ei a J\\ (26) 

for any function f(U). 



Note that the "coupled cluster equations" (23) and (24) are still rigorous. 

The linked cluster structure of the functions S(U) and F(U) follows from the fact that they 
may be expressed with the help of the projection operators (12) and (13) as in (15) in terms of 
"intrinsic" functions. These intrinsic functions are given by linked clusters and are defined as 
follows. 

Suppose that F^/q describes a state with Euclidean quantum numbers (p, T, v, v') (see eq. 
(15)). We then write F and S in the form (S has to have trivial quantum numbers) 

F(p,T,v,v') = Il lt ( P )U lr (T;u,u')F mt ( P ,T,u,u') (27) 

S = TloSint 

n = n zt (o)n, r (o,o,o) 

If F^/ corresponds in the continuum limit to a bound state, we expect that F int may 
be chosen to describe a localized state. This is analogous to non-relativistic many-body theory 
where bound states can be separated into square integrable functions of the relative coordinates 
and an overall center of mass motion, described here with projection operators. 

In analogy to nuclear matter for instance, the same localization holds true for the vacuum 
function Si nt because correlations have a finite range. 

The validity of these properties of the intrinsic functions is seen below through the structure 
of the expansion of these functions in terms of a localized basis, i.e. a basis of linked clusters. 

We shall first characterize this basis through its general properties and then describe the 
concrete construction in section 4. 

We call the basis 

X a iU h ,..,U lma ) a = 1,2,3,.. (28) 

and impose the following conditions: 

1) x a should be gauge invariant. 

2) x a should be "linked", see section 4 for the precise definition. A main consequence is 
that m a is finite for any a, though not limited. 

3) x a should be "standardized" , i.e. the equation 

T(u) X a = \X P {u e G E ) 

should only have solutions for a — (5. 

4) x a should be a strong coupling eigenfunction, i.e. 

E l,aFl,aX a = ta,lX a 

a 

It will be convenient to specify x l as the "plaquette function" by putting 

IIox 1 = 4(D - l)V 
and to distinguish the constant function via 

X° = l 

X° fulfils the relation 

n X° = s Nx° (33) 
where the symmetry factor so = \Gi r \ is equal to 8,48 for D = 2,3, respectively. 



(29) 



(30) 



(31) 



(32) 



Simplifying (27) by writing F = HF int and introducing the functions S(S) and F(S, F) by 



n s (34) 
(nF int )^(n s int )^ = iiF 

the coupled cluster equations (23) and (24) can be rewritten as 

{Smt) ^ + S ~ 4{D~T) X 1 = S °^k ( 35 ) 
{Fint)iin + 2F = (E — E )F int 

The "linked cluster theorem" for our lattice Yang-Mills theory consists then in the statement 
that if (for iV — > oo) we have the norm relation \S int \ : \F int \ oc 1, we have also \S\, \F\ oc 1. 

The norm relations for Si n t and Fi nt are fulfilled because of their localized nature, those for 
S, F follow then as a result of the fact that the "derivative" S^F^ selects only linked clusters. 
We shall prove this structure by suitable expansions in terms of the linked cluster basis (28). 

Introducing the summation conventions 

E = E (36) 

ct a=l,2,... 

E = E 

a a=0,l,2,... 

the expansions of the intrinsic wavefunctions read (specifying again the Euclidean quantum 
numbers) 

S in t(U) = J2 S «X a (U) (37) 

a 

F mt (p,T,u,u';U) = ]TF a (p,r>,z/V(f/) 

a 

The coupled cluster equations (23) and (24) may then equivalently be formulated as equa- 
tions for the coefficients S a and F a (p, T, v, v') 

e a F a (p,T,u, v') + 

2 ]T c7^(p,r,^y)^F 7 (p,r, u, Vl ) = (e- £ )F a (p,r>y) 



I 

The crucial quantities in these equations are the coupled cluster matrix elements Cf 7 and 
Cf 7 (j9, T, vi, u') defining the expansion of the functions S and F, respectively, which are ob- 
tained as follows: 

Determine first the set of numbers c^ 7 (u £ Ge) defined by 

E X^(THx 7 ), = E c ?« T Hx a (39) 



This yields then 



cf( P) r>y) = £ c^jUR- 1 ^- 1 ^ (40) 

u=(R, a) 

cf = Cf(0, 0,0,0) 

A proof hereof is given in the appendix. The linked cluster theorem guaranteeing the correct 
volume dependencies of the relevant quantities discussed above is now given by the fact that, 
due to the localized nature of the functions x a i t ne r.h.s. sums in eqs. (39) and (40) run only 
over a finite number of terms. 

The coupled cluster equation (38) also have the property that approximate solutions gen- 
erated by truncations (see e.g. Refs.[|, 0]) display correctly all relevant volume dependencies. 

The important task is now to compute the matrix elements c^, for which one needs an 
efficient and systematic way to set up and handle the basis elements \ a - 



4 Construction of the loop space basis 



In principle, an orthogonal basis of the type \ a has been constructed in Ref.[|ll[]: for simplicity, 
we will formulate the method for SU(2) and D = 2, but the generalization is obvious, though 
technically more difficult. 

Suppose x a ( a fixed) depends on the link variables (Ui, .., U r ). This basis function is then - 
- up to a normalization factor - uniquely characterized by the following set of angular momenta: 

1) We have a set (Ji, .., J r ), i.e. one (half integer) angular momentum for each link . 

2) We have an angular momentum J ab = J cd for each quadruplet of links l a ,--,ld forming 
a 4-point vertex in the link pattern (Z 1; ..,l r ). Here the convention is that the links (l a ,h) are 
oriented such that they are going into, and (l c , l^) are leaving, the common site. 

These angular momentum quantum numbers are constrained by J a = J& if (l a , lb) form a 2- 
point vertex and by the condition that the coupling J a + J& — »■ J c should be possible if (l a , If,, l c ) 
form a 3-point vertex. 

For instance, putting r = 7 and choosing the link pattern of Fig. 1, one has just three angular 
momenta (Ji(— J2 = J3), J&(= J5 = J&), J7) yielding the basis elements (a = (Ji, J4, J7)) 

X = E ^ 1 ,M 2 (^l)^ 2 ,M 3 (^)^ 3 ,M 8 (^ 1 ) (41) 
Mi,..,Mi 



Jl,J7,Ji \ D J 7 (TT ) ( ^ 1 '^7)«^4 

M 8 ,M 10 ,M 4 i M7 ' Ml ° l 7) \M 1 ,M 7 ,M 9i 



D A1 7 ,M 10 ( U 7) 

Dm 4j m 5 (U 4 1 )Dm 5:M( .(U 5 1 )Dm 6 ,m 9 {Ui 
We call this orthogonal set of functions the "D-loop basis" . 



In Ref. [12] this framework was used for estimating observables within the ELCE method, 



but higher order calculations were limited by the necessity to handle an increasing number of 



SU (n) couplings. Also it should be mentioned that in Refs. fTT , |I2f the linked cluster form of the 
lattice Yang-Mills many-body problem was taken into account within a different computational 
framework. 

An alternative for the construction of a basis is related to an exp S generalization of the 
Lanczos idea and was pursued in Refs. [6-10] for trivial representations of the Euclidean group. 
We will define this method here in such a way that it allows the computation of arbitrary 
Euclidean representations and also a transition to the independent, orthogonal D-loop basis. 
This yields especially a systematic way of eliminating linear dependencies. 



Starting with = x 1 from (31) we define "character functions" .., 4> s,ns by the iterative 
condition that the following expansion should hold (we put <ft 0,1 — l,no = n\ — 1) 

^(T(u)^ k % = £ £ C 5 sC 5 v k2T (v)^ k (42) 

v^G E 0<S<&i+5 2 ,k<n 5 

ki < n Sl k 2 < ns 2 

Because of the "derivative" /i, the functions <p 5,k are by definition linked - both terms on the 
l.h.s. have to have a common link variable for a non-vanishing result - and they also can be 
chosen to be standardized - i.e to obey the condition (29) - because we included the Euclidean 
operator T(v) on the r.h.s. of (42). 

The definition of the functions <fi s,k is made complete and unique by the condition that it 
should be just a product of characters, i.e. for each (S, k) there should exist a set of Li, .., L r 
defined by the generalized links 

Lj = (Iji, <Tji, .., lj mj , Oj m . ) (43) 

such that 

<y J - n; ,/r ( r;v. ./•;;;;;) (44) 

Hereby, for SU(3), the loops L±, .., L r should have all possible orientations compatible with the 
standardization of (p s ' k . For SU(2), however, all loops should have the same (fixed) orientation 
which is no loss of generality because of the relation trg = trg^ for g e SU(2). 

The result of the l.h.s. of (42) may be expanded in such terms because of an inductive 
argument: has the form (44). Assuming the form (44) for the two terms of the l.h.s. of 
(42) which we call, specifying for simplicity only the dependence on a certain (common) link 
variable Ui = V, tr(AV) for the first term and tr(BV a ) for the second term (a = ±1), the 
"differentiation" with respect to fi = (/, a) may be evaluated using eq. (19) and the standard 
property of the SU(n) generators A a : 

EW = 5(M^-;W (45) 

a Z n 

yielding 

Y,{tr{AV)) ha {tr{BV°)) ha = a{\tr{V»AV p B) - ^tr{AV)tr{BV a )) (46) 
a I In 

where p = (cr + l)/2. If the variable V = U\ also occurs in functions A{U) or B(U), additional 
terms arise on the r.h.s. due to the product rule of differentiation, but these terms will again 
display a loop space structure of the same type. The same happens if the first term has the 
form trlAV' 1 ). In this sense, (46) describes the "generic" case. 

Each character function <p s ' k , constructed in this way, is uniquely characterized by the set 
of "geometric" loops indicated in (44). For SU(2) and D = 2, examples up to third order are 
given in Figs. 2 and 4. For this dimension and for SU(2), the number of loop space functions 
is 4, 16 for the orders 5 = 2,3, respectively. 

5 The D-Mapping 

Of course, the system of loop space functions cf> s ' k is neither orthogonal nor linearly independent. 
For SU(2), for instance, the number of independent functions is known to be respectively 1,3,10 
up to third order (?J. 



However, the construction yields directly the expansion coefficients needed in (39) by taking 
in (42) the sum over u G Ge- 

The main problem which remains is to select an independent subset in the space of the 
functions (p s,k . 

In Refs. |7], [T(J independent functions were determined using besides trg = trg^ the relation 
tr(gg') = trg{trg') - tr(g^g') for g,g' G 377(2). 

The emerging functions were in general neither orthogonal nor unique. 

Within this paper, we propose a different strategy, namely, to relate the functions <fi s ' k - 
divided into convenient subsets - directly to the orthogonal "D-loop basis" \ a by a characteristic 
mapping, called D-mapping. 

This allows us to do the final calculation, i.e. to solve (approximately) eq. (38), with respect 
to the D-loop basis. But, at the same time, the crucial matrix elements c^ u may be computed 
in terms of the character functions using the D-mapping, while avoiding any explicit SU(n)- 
coupling or recoupling. We hope that this simpler structure will finally allow calculations of 
the type of Ref.[|7|, ^ [TTJ to higher order and/or for D = 3. 

The construction of the (non-invertible) D-mapping relies on the following structure of the 
D-loop basis already indicated in section 4. For any \ a — X a (Uh, ••, Ui ) (a fixed) there exists 
a maximal set of commuting operators (At, ..,A Ma ) with the property that 

A xX a = a xX a (47) 

and such that the state \ a is uniquely characterized by the eigenvalues (ai, .., omJ. 

In the case SU(2) and D = 2, these operators are apparently given by a combination of the 
two sets 

E\ a E\ a A = h, .., l ra (48) 



and all operators 



J2(E ba + E ca )(E ba + E ca ) (49) 



which fulfil the condition that (b, c) are outgoing links of a four point vertex in the link set 
related to x a - 

The more general case may be extracted from Refs. |i~T, 12|. Important for our purpose is that 



each A\ is a Casimir operator of the local lattice gauge group Gi oc , i.e. it is a certain polynomial 
in the operators E^ a of the type given above. For SU(3), two generalizations have to be taken 
into account: The third order Casimir operators have to be added to the set (47). In addition, 
suitable permutation operators have to be included if the SU(3) Clebsch Gordan decompositions 
generalizing eq. (41) have the property that the same irreducible representations occurs several 
times. (For SU(2), an example for the definition of such a permutation operator is given in the 
appendix A3.) 

The main point for the construction of the D-mapping is that the evaluation of the operators 
A\ on the states (p s,k can be done in precise analogy to the computation in eq. (46). The 
important ingredient is again the relation (45) and a corresponding third order generalization 
for SU(3) (see, e.g.|]i~7|). In other words, it is possible to compute the matrix elements of A\ 
defined by 

A x ^ k = £ {A x )f^ (50) 

8><S,k> 

These matrices related to A\ are finite (5' < 5) and stay small, in general. This is because the 
terms on the r.h.s. have to be consistent with the loop pattern of the variables of the l.h.s.. 
Since link variables may be removed by A\ (see eq. (46) for a = —1), this consistency also 



allows terms with removed variables on the r.h.s.. In the classification of the states as D-loop 
functions, this corresponds to the possibility that one of the link angular momenta may be zero. 
The loop pattern of (f) 5 ' k also determines the choice of the possible operators A x in (45). More 
details are given in section 6 where it is also shown that with the knowledge of the matrices of 
A\ it is sufficient to work out (42) without the derivatives. 

Eq. (50) corresponds to an evaluation of the operators A\ with respect to a non-orthogonal 
and overcomplete basis. Because the A\ are hermitian and commute with each other, they 
must nevertheless be simultaneously diagonalizable. Therefore, there must exist combinations 

v 7 = Eq/ fe (51) 

S,k 

such that 

A A y? 7 = a^y? 7 A = 71, .., 7 r7 (52) 

The choice of the operators A\ is determined by the loop pattern of the variables occurring in 
^. 

Any state ip 1 is by construction proportional to a D-loop basis function characterized by 
the eigenvalues a 7 ^, .., a 7 .. 

Consequently, the states •p> 1 are - up to a normalization factor and a possible Euclidean 
mapping (11) - equal if and only if their eigenvalue patterns are equal. 

Some subset of the functions y? 7 are then independent and orthogonal. They fulfill all 
conditions of (28) and can be identified with a certain subset of the x a - The D-mapping is 
just the restriction of (51) to such independent solutions and to the computation of the relative 
normalization factors for the dependent states y? 7 '. Having determined the D-mapping, eqs. 
(51) and (42) contain all ingredients for the computation of the crucial matrix elements (39). 



6 Computational strategy and examples 

We now describe the computational steps which - put into the language of a suitable computer 
program - would lead to the possibility to determine approximate glueball spectra. We shall 
elucidate these steps by some examples of low order for D = 2. 

1) Set up the character functions. 

In order to minimize the computational effort, we propose to divide the character functions 
into subsets of the following type. 

Introduce first the set of "generic" functions 

A^ fc 5 = 1,2,.. ; k = l,..,n s 

of 5-fold linked, standardized plaquette products. 

For SU(2), all plaquettes should have the same orientation. 

For D = 2, we have for SU(2) n s = 1, 2, 4, for SU(3) n s = 1, 4, 12 up to 5 = 3. Fig. 2 gives 
the corresponding loop patterns for SU(2). 

The relevance of this set of functions is two-fold: 

i) They determine the possible elements of the D-loop basis occurring up order 5. They are 
given by the link patterns of the generic set and the coupling rules of as many fundamental 
representations (and its adjoint) as there exist common links. Hereby, a double counting with 
lower order states has to be avoided. Fig. 3 exemplifies the related elements of the D-loop basis 
for SU(2) and D=2 up to 5 = 3. 



(53) 



ii) The set (53) is "generic" because each element defines a characteristic subspace given by 
functions 

A 5,k ' u ; v = l,..,M(6,k) (54) 

which is left invariant under the action of any Casimir operator of the lattice gauge group. 
Hereby, of course, only a finite number Casimir operators is relevant for any given (<5, k). The 
set of character functions (f) 5 ' K is contained in the set (54), (i.e. for a given 5 there exists for each 
k a pair (A;(/t), ^(/t)) such that that <p s > K = A S M^)A^ an explicit construction is not necessary 
because it is more convenient to work with (54). 

2) Compute the Casimir operator matrices. 

The set (54) is generated by applying the relevant Casimir operators on (53) yielding the matrix 
elements (50) as a system of - in general small - submatrices 

M(8,k) 

AxA s,k,» = J- A^(\,5,k)A 5 ' ky (55) 

v'=\ 

Applying the product rule and (46) for the evaluation of the Casimir operators yields for the 
loop structure of the subspaces (54) (5, k fixed) the simple geometrical condition that they are 
generated from the plaquette systems (54) by "cutting and glueing" doubly occuring links. Up 
to third order, the related SU(2) loop structures for D = 2 are given in Fig. 4, some examples 
of the corresponding Casimir operator matrices are presented in the Appendix A2. 

Note that the sets (54) may also contain elements of lower order if they occur during the 
cutting and glueing procedure. Also a standardization is not done. This is convenient since 
this makes the "Casimir matrices" (55) especially simple. 

3) Fix the D-Mapping. 

The next step is the diagonalization of the Casimir matrices (56) giving eigenfunctions 

^(6,k) = J2C:,(5,k)A s > k >»' (56) 

obeying 

A x ip v (8,k)=a x (v,8,k)<p v {8,k) (57) 

For the construction of the D-mapping one first may put with a suitable enumeration 

a(v, S, k) 

ip v {5, k) = N(u, 5, k)T(u(v, 5, k)) X * {vAk) ■ (58) 
Hereby, linear dependencies are eliminated by the identification prescription 

a(v, 5, k) — a(v', 5', k') a\(v, 8, k) — a\{v' , 5', k') for all X . (59) 

Of course, equality of the eigenvalue patterns guarantees the equality of the corresponding 
eigenfunctions only up to a (non-zero) factor N(u, 5, k) and up to a Euclidean transformation 
T{u{y,8,k)) 

For the computation of the normalization factors N(u, S, k) we observe that within our exp S 
framework (including a possible truncation) it is not necessary to work with basis states which 
are normalized to one. Hence only the relative factors are needed, i.e. we may put N(y, 8, k) — 1 
if the D-loop function x^' 6 '^ occurs for the first time when increasing the order 5. Also we 
may set for this first case u{v, 5, k) — 1. As a result we may find for each quantum number a 



an eigenfunction (56) characterized by (5(a), v(a), k(a)) defining an expansion of the elements 
of the D-loop basis in terms of the character functions 

x a = Y / C P J a \5(a),k(a))A S(a)Ma)y ■ (60) 

For each (5, k), the matrices C",(5, k) may be inverted, yielding with (58) the inverse mapping 
A*'*-" = Y / D u u ,(5,k)N(v',5,k)T(u(v',5,k)) X a(u ' Ak) ■ (61) 

Equation (61) - together with the inversion (60) - constitutes the D- mapping in a form which 
is sufficient for the computation of the coupled cluster matrix elements (39). 

We still have to give a recipe to compute the normalization factors N(u, 5, k). In principle, 
they could be determined by evaluating Haar measure integrals. This can be avoided, however, 
by rewriting the states <p u (6, k) in a (up to the normalization factors) unique form by using the 
usual procedure of eliminating linear dependences via the Cayley Hamilton relations 

trg^ = trg (62) 

g 2 = gtrg-1 (g e SU(2)) (63) 

and 

trg 2 = (trg) 2 - 2trg ] (64) 

g 3 = gHrg-gtrg^ + 1 (g G SU (3)) (65) 

This allows us to introduce a standardization of the functions A s,k ' u by eliminating for SU (2) 
(577(3)) all structures of the type trg n g' with n > 2 (n > 3). For this purpose, eq. (61) 
respectively (63) have to be iterated yielding formulas of the type g n = ag 2 + bg + c where a, b, c 
are polynomials in trg and trg^. (For SU(2), a = and b, c become polynomials in trg only.) 

For SU(3) and n = 2, also terms of the type trg 2 in (44) may be standardized with the help 
of (64). 

For the examples where the Casimir matrices (55) are computed, we give in the Appendix 
A2 also a construction of the corresponding part of the D-mapping. 

4) The incorporation of the Euclidean group. 

For the computation of the matrix elements (39) it is sufficient to work out (42) disregarding 
the derivatives, i.e. to determine the coefficients V^l^v ( u i v e Ge) given by 

A^T(u)A^ = J2vlCu% T (vW 3 (66) 

73,-" 

where we introduced the abbreviation (5, k, v) = 7. Here, the character functions A 71 and 
T(u)A. 12 have to fulfil the restriction that they are linked, i.e. they should have a common link 
variable. Also the trivial function A 0,1 ' 1 = x° should be left out. 

The following structures simplify the determination of these 77-coefficients: 

a) There is only one non- vanishing term on the r.h.s. of (66). If rp^^y is non- vanishing, it is 
equal to one. In this case we call the corresponding states A 71 , A 72 , A 73 non-trivially connected. 

b) For each triple A 71 , A 72 , A 73 of non-trivially connected character functions we have a 
characteristic set of Euclidean group elements u\,v\ such that 

A 71 T( MA )A 72 = T(v x )A 13 ; A = 1, .., n( 7 i, 72, 7s) (67) 



The determination of these elements U\, v\ is now simplified by the following structure: Sup- 
pose we have found all solutions u\,v\ for a non-trivially connected triple of generic functions 



Ag' fel 2> A )Ag' fe2 = 2> A )Ag' fcl (68) 

If v\ is suitably chosen, we have then for each A <5l ' felAl and A S2 ' k2 ' K2 a function A 5a ' k3 ' K3 so that 
they are non-trivially connected with the same set of Euclidean group elements as in (67) and 
this exhausts all possibilities. 

Given the generic functions of the r.h.s. of (68) and Ki,K2, the third character function 
cf) S 3' k 3' K 3 j s then determined by finding just one pair (uv) solving (65). 

Up to third order 5 3 = 3, a full computation of all //-coefficients in presented in Appendix 

A4. 



6) The computation of the c-coefficients in (39). 

Having solved the "combinatorial" problem of determining the coefficients of (66), one may 
compute the quantities c^\f* 2 by writing (66) in terms of the orthogonal and independent basis 
X° with the help of (60,61) and by applying the (Euclidean invariant) "total Casimir operator" 
J2ia EiaEia on both sides of the emerging equation. Writing 7(a) = (8(a), k(a), v(a)), we obtain 
as final result for the crucial coupled cluster matrix elements (39) 



c%:T = E(w + w-w)lE E ( 69 ) 

l u6G B 7li72,73 

°71 °72 iV W 4 ^ 73; ^(s( 74 ))-i^ 7 4 

74 with a (74) =03 



7 Discussion and conclusion 

The coupled cluster formulation of Hamiltonian lattice QCD needs an efficient method to 
deal with suitable basis systems of loop space functions. Within this paper we have demon- 
strated that it possible to combine the merits of a D-function basis, used within the ELCE 
framework |12| with those of the character sets used within recent coupled cluster attempts |7], || 



without facing the respective deficiencies. 

The merits are the orthogonality of the basis in the first case, the close relation to the 
Lanczos method and the easy computability of the coupled cluster matrix elements in the 
second case. 

The deficiencies are the need of handling too many SU(n) recoupling coefficients for the 
computation of the Hamiltonian matrix elements when using D-functions, the non-orthogonality 
and linear dependence of the states when using the character functions. 

Our combination is based upon the simple idea that the D-function basis may be charac- 
terized by the quantum numbers of a complete set of commuting operators. These operators 
are the Casimir operators of the local lattice gauge group (for the gauge group 577(3), also 
certain permutation operator have to be included) and our method relies on the fact that these 
commuting operators (where only a finite set is relevant for any specific case) may be evaluated 
as finite matrices with respect to the character functions. This allows the construction of a 
systematic mapping between the two frameworks. 

Invoking the lattice Euclidean symmetry of the regularized gauge field theory and system- 
atizing the action of this symmetry group, we were also able to formulate the coupled cluster 
lattice Hamiltonian eigenvalue problem for eigenstates with arbitrary lattice momentum and 
lattice angular momentum. The whole formulation may be done in the infinite volume limit. 



For any concrete calculation of the spectrum, a truncation prescription has to be defined. 



This point has been much in dispute [0, |7|, §, [L0| because previously one had to make a (non- 



unique) choice of independent functions from the non-orthogonal set of character functions. 
Within our method we have a more natural definition because the orthogonal D-function basis 
is uniquely determined. 

We want to stress that any truncated coupled cluster calculation will have the same limi- 
tations as any (finite volume) standard lattice Monte Carlo computation, namely that at best 
one has to hope for a scaling window indicating consistency with respect to the predicted renor- 
malization group structure which has to be displayed by any observable when approaching the 
continuum limit. (This structure is still unclear within Ref. flTol which gives the "best" coupled 
cluster results up to now.) 

The reason for this expected scaling window in given by the fact that the truncation which 
has to be defined with respect to an expansion of the intrinsic wave functions of the vacuum 
and of the hadron, necessarily limits the possible lattice volume over which the physical states 
may extend. Consequently, when the physical lattice scale is set by choosing the coupling g, 
the method has to break down when the physical lattice volume, given by the truncation - or 
by the number of lattice points in the standard lattice Monte Carlo case - becomes smaller than 
the size of the hadron. 

First attempts at doing concrete numerical calculations within the reported framework are 
on the way and will be reported in the future |18| . 



Finally we want to mention that our computational framework may, in principle, be easily 
extended to include Fermions. Especially, a formulation for Wilson Fermions within a quenched 
approximation yields equations whose treatment appear to be no more complicated than that 
for glueballs. Details of this structure will be reported elsewhere. 
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Appendix 

Al. Proof of equation (37) 

Introducing the abbreviation 

D^(R } a)=dlARy^ (70) 
the expansion of F^S^ according to (27) and (34) yields the relevant terms 

{Ii u { P )Ii lr {Y-^v') X a )^X P ), (71) 
= E D^iuJiTiujx-UTiuJx* 3 ), 

= E D ly(ux)T(ui) E X a »(T{u)x p )»\ 

U\ u 

= E^jK) T K)EC T Wx 7 



= E E («- 1 )T( UlU ) X ' 5r 

U,7 Ml 

= E HttdOiMr; ^i) E ^(« _ V 

i/l n,7 

Here we assumed (39) and we arrived at the result (40). 
A2. Examples for the D-mapping 

Within this appendix we will give the construction of "Casimir matrices" (55) for SU (2) and 
D = 2 for some typical cases of subspaces (5, k) taken from in Fig. 4. Subsequently, we will 
present the corresponding part of the D-mapping resulting from a diagonalization. 

Note that also the D-loop basis is characterized by the same quantum numbers (5, k) , see 
Fig. 3. 



(S,k) = (2,1) 

Here, two parallel links with the same orientation occur. If a given link of this type is denoted 
by U, this case yields for the corresponding Casimir operator (48) a 2 x 2 matrix with repect 
to states of the type 

A 1 = tr(BU)tr(CU) (72) 
A 2 = tr(BUCU) 

For (5, k) = (2, 1), we have B = C and U may be any of the four links. A more general case 
is e.g. given by U — U2 for (S,k) = (3,2). (For enumeration see Fig 3). The evaluation of 
A = E a E a yields 

A A 1 = A x + A 2 (73) 
AA 2 = A x + A 2 

The diagonalization gives the Casimir spectrum as an example of eqs. (56) and (57) 

cp 1 = A 1 + A 2 ; a(l) = 2 (74) 
if 2 = A 1 -A 2 ; a(l) = 

For A = B one may eliminate in the state f 2 the function (AU) 2 via (61) yielding the identifi- 
cations with the D-loop functions (see Fig. 3 for notations) 

y? 1 = x 2 ' 1 (75) 
V 2 = 2 = 2 X ° 

Note that we put the normalization factor N(u, 5, k) equal to one when the corresponding basis 
state x a occurs for the first time. Within our examples, the Euclidean mapping T(u) in (58) is 
the identity in most cases. 
With the inversion 

A 1 = Ix^ + X (76) 
A 2 = l^-X 



these formulas give the D-mapping relevant for the subspace (5, k) = (2, 1). 



(<5,k) = (2,2) 

This case is of the type 



yielding 



A 1 = tr{BU)tr{U ] C) 

A 2 = tr(BC) 

AA 1 = 2A 1 -A 2 

AA 2 = 



A diagonalization of (80) gives 

V? 1 = 2A X -A 2 ; a(l) =2 



ip 2 = A 2 ; a(l) = 



With the identification 



and the inversion 



v? 2 = x 2 ' 2 ' 2 



A 1 = ^(f-' + f 2 ) 



A 2 = X 2 ' 2 ' 2 

this defines all the D-mapping ingredients for (5, k) = (2, 2). 
(5,k) = (3,1) 

Here one has to deal with a three-dimensional subpace given by 

A 1 = itrgf 
A 2 = trg trg 2 
A 3 = trg 3 

with g = BU . The product rule yields for the evaluation of the Casimir operato 

A A 1 = -A X + 3A 2 
4 

AA 2 = A X + -A 2 + 2A 3 
4 

AA 3 = 3A 2 + -A 3 
4 



The eigenvectors and eigenvalues are 



15 
T 



V? 1 = A X + 3A 2 + 2A 3 ; a(l) 
V? 2 = A 1 - A 3 ; a(2) = | 
ip 3 = A 1 - 3A 2 + 2A 3 ; a(3) = -jj 



Using trg 3 = (trg) 3 — 3trg obtained from (63), this yields the identifications 

V 1 = X 3 ' 1 (85) 

Since the eigenvalue a(3) is negative, we must have 

V? 3 = (86) 

yielding the linear dependence relation 

A 2 = i(A 1 + 2A 3 ) = (trgf-2trg (87) 
3 

which is just the result for A 2 when eliminating g 2 by (61). The inversion (64) now reads 

A 1 = ix 3,1 + 2x M (88) 

A 2 = U"' 1 
6 

A 3 = ^-X 1 ' 1 



(5,k) = (3,4) 

This case is interesting because it involves a 4-point vertex, the dimension of the subspace (54) 
is 5. Calling the doubly occuring link variables U~i, U 2 the generating character states are of the 
type 

A 1 = tr{BU{} tr{U\CU 2 ) tr{U ] 2 D) (89) 

A 2 = tr(BCU 2 ) tr(U ] 2 D) 

A 3 = tr(BUi) tr{U\CD) 

A 4 = tr(BCD) 

A 5 = tr{U\CU 2 ) triBU^D) 

Now we have three relevant Casimir operators, A\ and A 2 as before of type (48) and A% = 
(E la + E 2a )(E la + E 2a ) - A 1 - A 2 = 2E la E 2a which is of the type (49). The related Casimir 
matrices are given by 

(90) 



A 1 A 1 


= 2A l 


-A 2 


A ± A 2 


= 




AiA 3 


= 2A 3 


-A 4 


A ± A 4 


= 




A ± A 5 


= -A 4 


+ 2A 5 


A 2 A l 


= 2A 1 


-A 3 


A 2 A 2 


= 2A 2 


-A 4 


A 2 A 3 


= 




A 2 A 4 


= 




A 2 A 5 


= -A 4 


+ 2A 5 



(91) 



A3A 1 = -2A 1 + A 2 + A 3 - A 5 (92) 



A 3 A 2 = 

A 3 A 3 = 

A 3 A 4 = 

A 3 A 5 = 2A 4 -4A 5 



The simultaneous eigenfunctions are 



V? 3 = -2A 1 + A 2 + A 3 -A 4 + A 5 ; ai(l) = a 2 (l) = 2 a 3 (l) = -2 
V? 4 = -2A 5 + A 4 ; ai(2) = a 2 (2) = 2 a 3 (2) = -4 
V? 2 = -2A 2 + A 4 ; ai(3) = a 2 (3) = 2 a 3 (3) = 
if 5 = -2A 3 + A 4 ; ai(4) = 2 a 2 (4) = a 3 (4) =0 

V? 1 = A 4 ; ai(5) = a 2 (5) = a 3 (5) = (93) 

with the identifications 

^ = x 3 ' 4 '" v = 1,2,3,4 (94) 

^5 = T(P)x 3,4,2 

where P describes the reflection (parity transformation) defined in equation (96). 
The inversion (64) reads 



A 1 = 




A 2 = 


l(x 3A1 -x 3A ' 2 ) 


A 3 = 


-\(T(P)X 3A2 - X 3A1 ) 


A 4 = 


x 3,4,l 


A 5 = 


\{x 3Al -x 3AA ) 



,3,4,2 ,,3,4,4\ 



(95) 



A3. The "local action" of the permutation group 

We explain this structure for the "typical" example (5, k) = (3, 4) discussed in the Appendix 
A2, a generalization for general cases is straightforward, but will not be displayed within this 
paper. 

For SU(2), first the equivalence of the fundamental representation and its adjoint has to be 
invoced by introducing the skewsymmetric 2x2 matrix 

£ = (-°i I) 

which has the property that eg^e -1 = g for any g G SU{2). Consequenty, the modified link 
variable U = Ue obeys instead of (5) the "tensor product" transformation rule (written in 
terms of matrix elements) 

(Ui)jk = g(x) jf g(x + ej) kk >{Ui)j> v (96) 



Introducing the "modified loop group elements" specified according to the common (four point) 
lattice site of our example (see eq.(89)) 

a = BUie (97) 
= U\DU 2 e 
7 = U\De 



the character states (89) may be written as 

A- = d(uy^^a n32 p 33j4lj536 ( 98) 

Each d-coefficient couples the tensor product of the six fundamental representations - - defined 
"locally" corresponding to the chosen the four point vertex - to the trivial representation, i.e., 
introducing the related six angular momentum operators s(r), r = 1, .., 6 (sj are the Pauli spin 
matrices in our case) the states (98) obey 

StotojA" = (99) 

6 

Stotal = S ( r ) 
r=l 

Obviously, the "total angular momentum" s tota i is invariant with respect to any permutation 
of the six variables appearing as indices in (98), i.e. we may simultaneously characterize the 
space (89) the representations of the permutation group S 6 which acts on the states (98) by 

^ h ,yii2hjiji36 _^ ^^ z/ y CT (i)i t T(2)i t T(3)i<T(4)i CT (5)i t T(6) 

(o" G Sq). One may use the decomposition of this representation for classifying the states 
in the space (89). In its general form, however, the corresponding permutation operators do 
not commute with the Casimir operators (47) since they involve only "reduced total angular 
momenta". In our case we have e.g. 

A 1 = (s(2) + s(3)) 2 (101) 
A 2 = (s(4) + s(5)) 2 

A 3 = (s(2) + s(3) + s(4) + s(5)) 2 -A 1 -A 2 

Conveniently chosen subgroups of Sq, however, do commute. We may take, e.g., S 2 embed- 
ded in Sq in different ways: If 7r is the non-trivial element of S 2 , we may put 

7r(l, 2, 3, 4, 5, 6) = (1,3,2,4,5,6) - yielding x 3 ' 4,1 an d X 3 ' 4 ' 2 as antisymmetric and x 3 ' 4,3 and 
X 3 ' 4,4 as symmetric representations - 

or 7r(l, 2, 3, 4, 5, 6) = (1,3,2,5,4,6) - yielding x 3,4 ' 1 an d x 3,4 ' 3 as antisymmetric and x 3,4 ' 2 
and x 3 ' 4,4 as symmetric representations. 

Of course, for SU(2), this does not yield independent quantum numbers. With a suitable 
choice of the permutation subgroup, however, this may be the case for SU(3). 

A4. Examples for the incorporation of the Euclidean group 

We restrict ourselves to SU (2) and D = 2. A convenient enumeration of the Euclidean group 
for D = 2 is given by 

[n, //, m, a] = R n t^R m P a ; m,n = 0,1,2,3 (102) 

a = 0,1 

pt = 0,1,2,3,4,5, 



(100) 



where the use the following conventions: 

The parity transformation P is fixed by the condition T(P)A 3 ' 4 ' 2 = A 3 ' 4 ' 3 

The rotation R is given by the constraint that it has rotation angle ir/2 and fulfils T(R)x l = 

X 1 = A 1 ' 1 ' 1 

The "one-unit" translation t is defined by the condition that Aq 2 = x l Tit^x 1 - 
Taking into account T(P)x 1 = X 1 ) we obtain in lowest order 71 = 72 = (1, 1, 1) the non- 
trivially connected cases (n, m, a are arbitrary with the restriction (96)) 

73 = (2, 1,1) u= [n,0,m,a] v = [0,0,0,0] (103) 
73 = (2,2,1) u = [ra, l,m, a] v = [4- n, 0,0,0] 

Combining first and second order on the r.h.s. of equation (69) we have the following generic 
alternatives: 

7l = (2, 1,1); 72 =(1,1,1): 

73 — (3, 1, 1) u = [n, 0, m, a] v = [0, 0, 0, 0] (104) 
73 = (3, 2, 1) u = [n, l,m, <t] v = [4 -n, 0,0,0] 



7l = (2,2,1); 72 = (1,1,1): 



73 = 


(3,2,1) 


u = 


[0, 0, m, a] 


v = 


[0,0,0,0] 






u = 


[0, 1, m, a] 


v = 


[0,1,2,0] 


73 = 


(3,3,1) 


u = 


[0, 2, m, a] 


V = 


[0,0,0,0] 






u = 


[2, 0, m, a] 


V = 


[0,1,2,0] 


73 = 


(3,4,1) 


u = 


[l,l,m,<j] 


V = 


[0,0,0,0] 






u = 


[3, 1, m, a] 


V = 


[0,0,1,1] 



Here, the Euclidean elements v are chosen such that a non-trivial connection with the same 
pairs (w, v) is described for the 7-tripletts 

(2, 1,2), (1,1,1), [(3, 1,2) or (3,2,2)] for the two cases (103) and 

(2,2,2), (1,1,1), [(3,3,2) or (3,4,2)] for the cases (104). 

The Euclidean elements (w, v) for the cases where 71 and 72 are exchanged may be obtained 
by the replacements 
u — > u -1 , v — > u~ x v 

in the above formulas. This follows from A 2 T(«- 1 )A 1 = T(«- 1 )(A 1 T(m)A 2 ) = T^-^A 3 if 
A 1 T(m)A 2 = 2»A 3 . 

This provides all non-vanishing 77-coefficients (66) up to the order 5 3 = 3. 
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Figure 1. Link pattern of the D-loop functions (41). The numbers indicate the enumeration 
of the link angular momenta (Ji, .., J 7 ). 
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Figure 2. Loop structure of the generic character functions A^ fc (see equation (53)) up to order 
5 = 3 for SU(2). 
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Figure 3. SU(2) D-loop basis functions x S ' k,u characterized by the link patterns and and 
the related possible Casimir eigenvalue patterns which are generated up the order S = 3. Non- 
equal angular momenta are enumerated in the link patterns and indicated in that order in the 
eigenvalue patterns. Upper indices stand for the degeneracy of these angular momenta. 

For 5 = 3, k = 3 or 4 all angular momenta which are not specified are equal to \. 

For (5, k) = (3,4), the third angular momentum is given by the "intermediate" coupling 
J1+J2 (eq. (49)). 

The last column gives the linear dependent D-loop functions emerging up to this order. 
The orientation of the links is not marked, it may be taken analogously to Fig. 1. 
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Figure 4. Loop structure of the functions A 5 ' k ' u (see equation (54)) up to order 5 = 3. 



